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| Abstract The Ricci curvature criterion is used for the investigation 

Q ■ of the relative instability of several configurations of N-body gravitating 

systems. It is shown, that the systems with double massive centers are more 
unstable than the homogeneous systems and those with one massive center. 
In general this shows the efficiency of the Ricci curvature method introduced 
by Gurzadyan and Kocharyan (1987) for the study of N-body systems via 
relatively simple calculations, i.e. for small N, and hence small computer 
resources. 



> 
(N 



O 

^ : 1 Introduction 



Oh: 

6 

H 



The interest to the problem of dynamics of stellar systems with double 
massive centre is due to the discoveries of double nuclei in the centers of 
several galaxies have various separation in projection - from 2 pc for M31 
"^1 up to 800pc for Markarian 273; double massive object was descovered in the 

centre of Arp 220 g] too. A massive binary system located in the core of a 
galaxy can lead to a number of dynamical effects. 
^ I In the present paper we investigate the role of double centers, i.e. binary 

$— i ' massive bodies situated in the center of N-body systems, in the dynamical 

instability of the gravitating systems using the Ricci criterion (has been 
introduced in M) for estimation of the relative instability (chaos) of those 
systems. This criterion, for example, had enabled to establish in that 
the regular central field is increasing the instability of the N-body systems. 
It is remarkable that this result has been obtained via numerical study 
of relatively small number of particles and later has been confirmed with 
extensive simulations on powerful computers [Q]. The effect of the central 
regular field is crucial for the understanding of the relaxation, mixing and 
evolution of the galactic cores. The situation is complicated in a sense that 
there are both type of effects - acting to increase and decrease the chaos in 
the system. 
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The flow of numerical methods range from the Lyapunov numbers and 
KS-entropy, approximate expansion and frequency maps up to powerful 
methods based on direct solutions ||, |J, |7j. The adequacy and efficiency 
of a given method in each particular case is itself an interesting problem. 
For example, the efficiency of estimation of Lyapunov numbers is limited 
by a number of reasons, particularly due to the exponential growth of the 
errors at large enough N inevitable at any iterated numerical procedure. 

Geometrical methods based on the theorems of the theory of dynamical 
systems provide an alternative way of study of N-body systems via reduc- 
ing the problem to that of the geometry of the phase and configurational 
spaces of the system ||. In physical problems this method has been initially 
used by Krylov || and for N-body gravitational systems - by Gurzadyan 
and Savvidy [10|. In the latter papers it was shown that spherical N-body 
systems are K-systems, i.e. are mixing systems with exponential instability. 

The statistical properies plays a key role for the understanding of the 
relaxation and evolution of many astrophysical objects - from the Solar 
system to galaxies and clusters of galaxies. 

First, we will briefly describe the Ricci curvature formalism, and the 
algorithm of the numerical calculations. 



2 The Ricci curvature criterion 

N-body gravitating system with potential 

N 

V = — G Grrij mj / rjj , 

i<j 

(rij is the distance between the particles with masses and mj) via a 
variational principle can be transformed to a geodesic flow in a Riemannian 
space@. The behavior of close geodesies in this space is described by the 
Jacobi equation 

V u V u n + Riem(n, u)u = 0, 

where u is the velocity of geodesies, n is the separation vector of close 
geodesies and V denotes the covariant derivative. For the normal com- 
ponent of the deviation one can obtain the following equation 




= -2K u . n || n || 2 +2 || V«n 
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where 

[Riem(n, u)u]n 

II n II 

is the so-called two-dimensional curvature (Riem is the Riemannian curva- 
ture). 

If K u ^ n is strongly negative in all two-directions (u,v) and everywhere 
in a compact manifold, the system possesses maximally strong instability 
properties is isomorphic to Bernoulli shift and is an Anosov system 

However this condition is too strong and therefore is not fulfilled for 
real physical systems. It is reasonable therefore to look for a weaker crite- 
rion using some average deviation of geodesies and a mean curvature in the 
manifold M. Consider the Ricci curvature 

u^u u 

II U II 

where is the Ricci tensor. 

The criterion of relative instability based on the Ricci curvature reads 
(Gurzadyan and Kocharian ||): 

The geodesies 71 (s) with velocity u\ is more unstable with respect to the 
geodesies 72(5) with velocity 112 within some interval [0; Si], if 

r = inf[r n (s)l 

m I uV )\ 

and 

T\ < r% T\ < 0. 

The advantage of this criterion is that it is checkable via computer simu- 
lations of many dimensional systems, including of N-body systems. Note 
also, that as distinct from the Lyapunov numbers, this criterion describes 
the local (in time) properties of the system. 

3 Numerical simulations 



At our computer experiments we estimated the Ricci curvature for several 
evolving configurations, i.e. traced the variation of the Ricci curvature in 
time. The formula for the Ricci curvature for N-body systems can be derived 
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from the formulae given in previous section and has the following form [0], 



(3N - 2) Wik-uW + 3 _ 2) _(WX) 2 (3JV - 4) \VW\ 



w 2 VF 4 ; 4 

where 

dW 

W = E-V: Wi 



dqi 



d 2 W 



W ik = ^——; AW = Y,Wu; 
oqiOq k ^ 

We built the systems with N = 22 using the scheme described in 0, i.e. 
we considered systems, so that the particles were located at the apexes of 
concentric cubes of unit sides. The velocities of particles were chosen in way 
to have no rotational momentum for the system. 

We estimated the variation of the Ricci curvature by time for the follow- 
ing configurations (Figures 1-3): 

1. Homogeneous, i.e. all particles have the same mass m; 

2. With one central massive particle mi, while N — 1 particles have the 
same mass m « M; 

3. With two massive particles of masses mi and m<i situated in the 
central part of the system, while N — 2 particles have the same mass m « 
mi,m 2 ;. 

Our calculations showed that the systems with massive center are more 
unstable than the homogeneous ones, thus confirming the conclusions in 0], 



H], Q, [14|. Note, the growth of the instability with the increase of the 
central mass M. 

Figures 1-3 illustrate the comparative instability of the three different 
types of systems mentioned above. The most unstable among the initial 
configurations is the system with double-central masses. Note, that as the 
systems evolves the Ricci curvature tends to zero for all three systems, how- 
ever the rate of tending is most rapid for the system with double centers. 
Physically it is clear that the third type of systems has to dissolve quicker for 
small number of particles. Just this tendency has been also noticed at the 
numerical experiments, namely the tending to zero of the Ricci curvature 
becomes slower with the increase of N. 
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In other words, the double massive central objects make the system 
more unstable initially, however then the system evolves quicker to its final 
dissolved state with regular orbits. 

This tendencies have been confirmed in numerous experiments with con- 
figurations with various initial conditions. 

4 Conclusions 

Thus we used the Ricci curvature criterion to study the relative instability of 
three types of N-body systems: homogeneous systems, those with one central 
mass and two central massive bodies. The following main conclusions had 
been drawn via the numerical experiments: 

1. The presence of the second massive central object makes the system 
more unstable as compared with that those of a single massive center and 
the homogeneous ones. 

2. The system with double massive objects is evolving more quicker 
towards dissolution, i.e. to a more global regular situation. 

3. The greater is the ratio of the mass of the massive objects to the mass 
of the rest particles, the more is the difference in the initial instability and 
the rate of evolution. 

The main conclusion is however the efficiency of the Ricci curvature 
criterion for the study of such complex many dimensional systems - N-body 
systems - by means of simple numerical experiments. 

We thank V.G.Gurzadyan and S. J.Aarseth for valuable discussions. This 
work is supported in part by INTAS grant. 
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Figure 1: Comparative instability of three types of systems: with equal 
masses of particles, with a single central mass ra\ = 10 and with two massive 
centres m\ = m,2 = 10; here and in the following figures mi, 777.2 are the 
masses of the central particle/particles, while the other TV — 1 particles have 
777, = 1. 
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Figure 2: Comparison of configurations with various parameters of masses. 
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Figure 3: Comparative instability of systems with higher values of central 
masses. 
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